frustration-induced phase transitions in the spin-S' orthogonal-dimer chain 
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We investigate quantum phase transitions in the frustrated orthogonal-dimer chain with an arbi- 
trary spin S > 1/2. When the ratio of the competing exchange couplings is varied, first-order phase 
transitions occur 2S times among distinct spin-gap phases. The introduction of single-ion anisotropy 
further enriches the phase diagram. The phase transitions described by the present model possess 
most of the essential properties inherent in frustrated quantum spin systems. 



Geometrical frustration in strongly correlated electron 
systems has attracted much current interest. A remark- 
able example of materials is the two-dimensional spin- 
gap compound SrCu2(BO3)2,0 in which the character- 
istic orthogonal-dimer structure of the Cu 2+ ions sta- 
bilizes the spin-singlet ground state. [J[ [|, [| Remark- 
ably, strong frustration induces an anomalous first-order 
phase transition in addition to the plateau-formation 
in the magnetization process. ]5| More recently, another 
orthogonal-dimer compound Nd2BaZn05^] was synthe- 
sized, where higher-spin and orbital moments (J = 9/2) 
show the antiferromagnetic order at the critical temper- 
ature Tn = 2.4K, making such higher-spin systems more 
interesting. 

One of the most prototypical phenomena inherent in 
the frustrated systems is the first-order transition, which 
is triggered by the competition of various phases due to 
strong frustration. Although the systematic treatment 
of first-order transitions is not easy, it is highly desirable 
to clarify how such first-order transitions are induced by 
frustration in order to understand the essential properties 
common to such frustrated quantum spin systems. 

In this paper, we investigate a remarkable one- 
dimensional (ID) spin-S 1 orthogonal-dimer model, which 
possesses most of the essential properties of first-order 
phase transitions in this class of frustrated spin systems. 
0, § ! By exploiting the non-linear sigma model 
(NLerM) approach as well as the exact diagonalization 
and the series expansion, we find the (25 + 1) distinct 
spin-gap phases in the spin-S* chain, which are separated 
by first-order quantum phase transitions. The effect of 
single-ion anisotropy is also addressed. We demonstrate 
that a higher-spin generalization of the model results in 
the remarkably rich phase diagram, which realizes the 
idea of valence-bond-solid (VBS) [JllJ in a sequence of 
first-order phase transitions. 

Let us consider the ID quantum spin system with the 
orthogonal-dimer structure || ||, |l0) shown schematically 
in Fig. [l]. The corresponding Hamiltonian reads 

H = J ^ (S2,j ■ Szy + S 3j ■ Sij+i) 

3 

+J'^(S 2j +S 4j ).(S 1j+ S 3j ), (1) 
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FIG. 1: Orthogonal-dimer spin chain: the dimer bonds shown 
by solid lines have the characteristic orthogonal structure. 



where Sjj is the i-th spin operator in the j-th plaque- 
tte and J(J') is the antiferromagnetic exchange coupling. 
The S = 1/2 model was so far studied well.|3, 10 Below, 
we focus on a higher-spin generalization of the system. 

We first exploit the NLerM technique to clarify the 
topological nature of the system. Jl^, [l3|, |lj, [ll], [l(| 
Note that the Hamiltonian has the remarkable rela- 
tion [H, S%j + &4:,j] = 0. Therefore, by introducing 



the composite spin Tj defined as Tj = S2.J + S^j, ' 
obtain the effective mixed-spin Hamiltonian as, H 
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1) — jJS(S + 1)N, where N is the total number of sites. 
Then, the Hilbert space of the Hamiltonian (Q) can be 
classified into each sub-space specified by [S; {Tj}]. The 
singlet ground state as well as relevant low-energy ex- 
citations are in the space with uniform Tj (= T). In 
particular, for a given T ^ we can describe low-energy 
properties by the NLctM (T = gives a trivial system 
with decoupled dimers). Introducing three kinds of the 
fluctuation fields, @ we obtain the Euclidean La- 
grangian C with the effective field <fi as, 
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where 6 = 2irT, g = 2AB/T, v s = QJ'StB/A with A = 
(2t + j ) 1/2 B = [4 J0 i 2 + 2(1 - 2j )t + 3o]- 1/2 , jo = 
J 1 1 J and t — S/T. Note that the topological angle 8 
is zero (mod 27r) for any choice of T, and the system is 
always gapped. Therefore, in case the phase transition 
between the distinct sub-spaces [S; {Tj}} occurs, it should 
be accompanied by the discontinuity in the parameters g 
and v s . This suggests that the possible quantum phase 
transition should be of first-order. 

To clarify this point, we numerically diagonalize the 
Hamiltonian ([!]) for a small cluster with periodic bound- 
ary conditions. The ground state energy is shown in Fig. 
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p|. It is found that the cusps appear 2S times in the 
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FIG. 2: Ground-state energy as a function of J'/ J obtained 
by the exact diagonalization for 8 sites. From up to down 
S = 1/2,1,3/2,2,5/2,3. The broken line is the energy for 
the classical limit S 1 ^ oo. The inset shows the derivative of 
the ground state energy, in which the origin for each curve is 
shifted for convenience. 

energy diagram for the spin-5 case, implying that the 
first-order quantum phase transitions indeed occur 25 
times. Also shown is the energy in the classical limit 
S -> oo, which is given as E/JS 2 N = -|[1 + (J' /J) 2 }, 
and — J' I J for the helical phase (0 < J'/ J < 1) and 
the antiferromagnetically ordered phase (1 < J' /J). As 
S increases, the profile of the energy diagram gradually 
approaches the classical one, although the 2S transitions 
should exist for any finite S. By examining several sys- 
tems with different cluster size, we determine the phase 
diagram rather precisely, as shown in Fig. pL In this 
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FIG. 3: Ground-state phase diagram for the orthogonal-dimer 
spin chain with generic spin S. The closed circles, open circles 
and open squares indicate the phase boundaries determined 
by the spin chain (N = 8, 16, 24) with periodic boundary 
conditions. 

way, the present model is the remarkable spin-5 micro- 
scopic model, which clearly determines a sequence of the 
phase transitions triggered by frustration. Moreover, the 
first-order transitions found here are contrasted to the 
second-order ones known for the ordinary spin-5 Heisen- 
berg chain with bond alternation. p3[ 



We now clarify the nature of each spin-gap phase by 
taking the S = 1 model as an example, which contains 
three distinct spin-gap phases. For this purpose, the VBS 
description |ll| of the ground state is useful, where the 
topological nature is specified by a combination of singlet 
bonds between the decomposed 5=1/2 spins, as shown 
in Fig. ||. Recall that in the small (large) J'/ J region, 
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FIG. 4: VBS picture of the phases for the S = 1 model. 

the composite spin T = 0(2) on each diagonal bond is re- 
alized. Hence the dimer (plaquette) phase characterized 
by Fig. U (a) [(c)] is stabilized there. On the other hand, 
in the intermediate phase, strong frustration induces the 
spin T = 1 on each diagonal bond, resulting in the singlet 
phase characterized by Fig. ^(b). This may be regarded 
as the frustration-induced Haldane phase. 

To confirm the above predictions based on the VBS 
analysis, we perform the numerical diagonalization of 
small clusters in the corresponding sub-spaces. In Fig. 
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FIG. 5: The energy for the S ■■ 
in the restricted sub-space. 



1 orthogonal-dimer spin chain 



H (a), the flat line is the energy for the exact dimer state. 
The energy for the Haldane phase (N = 8,16 and 24) 
obtained by the exact diagonalization in the sub-space 
\{Tj = 1}] is shown as the solid line with closed cir- 
cles, open circles and open squares, respectively. We also 
show the ground state energy of the original orthogonal- 
dimer spin chain (TV = 16) as the bold line with closed 
triangles. As seen clearly, the increase of J'/ J triggers 
the first-order quantum phase transition from the dimer 
phase (Tj = 0) to the Haldane phase (Tj — 1) at the 
critical value (J'/J) c ~ 0.56. Shown in Fig. || (b) is 
the energy around the other first-order critical point be- 
tween the Haldane phase (Tj = 1) and the plaquette 
phase (Tj = 2). In the plaquette phase, we have used the 
series expansion JlTj by choosing an isolated plaquette 
as the unperturbed system and regarding the interaction 
between plaquettes as a perturbation. The ground state 
energy calculated up to the ninth order is shown as the 
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broken line in Fig. from which we find the critical point 
(J'/J) c ~ 0.88 between the Haldane and the plaquette 
phases. It is thus concluded that the first-order quantum 
phase transitions occur among three singlet phases spec- 
ified by the distinct sub-spaces [S; {Tj}](Tj = 0, 1, 2), as 
predicted by the NLerM approach. 

Keeping this in mind, we now consider excitations in 
the 5" = 1 orthogonal-dimer spin chain. In the dimer 




FIG. 6: Various excitations in the S — 1 orthogonal-dimer 
spin chain. Fhe bold lines indicate a dispersive magnetic ex- 
citation, while the solid (broken) lines dispersionless magnetic 
(non-magnetic) excitations. 

phase (J'/J < 0.56), the ground state is given by the 
product of isolated dimers [{Tj = 0}], which allows us 
to estimate the excitation gap exactly from the finite- 
size calculation. The lowest magnetic excitation is de- 
scribed by a defect in the uniform spin alignment, i.e., 
it is given by the lowest triplet state in the space of 
[1, 0, 0, 0, ■ • ■]. It is also found that in the vicinity of the 
critical point (J'/J ~ 0.56), a non-magnetic excitation 
belonging to the space [1, 1, 0, 0, 0, • • •] can be the lowest 
one, as shown by the broken line. In the Haldane phase 
(0.56 < J'/J < 0.88), we estimate several kinds of spin 
gaps by extrapolating the results for N = 16, 24 in the 
formula A(N) = A(oo) + aN~ 2 . The Haldane gap ex- 
pected naively is the lowest away from the critical points, 
shown as the bold line in Fig. |6[ There are other disper- 
sionless excitations, which can be described by a defect 
in the spin alignment such as [0, 1, 1, 1, • • •], [2, 1, 1, 1, • • •]. 
These excitations are bound into another non-magnetic 
state, which is the lowest excitation in the Haldane phase 
close to the plaquette phase (J'/J ~ 0.88). In the pla- 
quette phase (0.88 < J'/J), the series expansion is more 
efficient to obtain the dispersive and the dispersionless 
excitations. The results computed up to the seventh and 
the ninth order are shown as the closed circles in Fig. ||. 
In this way, several different excitations become almost 
degenerate around the first-order phase transition points, 
reflecting strong frustration. 

Finally we discuss the effect of single-ion anisotropy 
with the Hamiltonian H D = DJ2 (Sfj) . The compos- 
ite spin Tj is not a good quantum number in the pres- 
ence of the D-term. The frustration-induced interme- 



diate phases, however, are stable against the introduc- 
tion of small anisotropy. In fact, the cusps still exist in 
the energy diagram for the spin S = 1 chain (N — 16), 
from which we determine the phase boundaries shown 
as the open circles in Fig. ffl. The Haldane phase grad- 
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FIG. 7: Phase diagram for the anisotropic S = 1 chain with 
the orthogonal-dimer structure. 

ually shrinks with the increase of D(> 0), and finally 
disappears. We note that around D > 0.8 the dimer 
and the plaquette phases indeed merge into the single 
phase. In order to clearly distinguish the Haldane phase 
and the dimer (or plaquette) phase, we make use of the 
symmetry of the space inversion P and the spin reversal 
t. Under twisted boundary conditions |fl8| (S® jy/4+i = 

- S i,f S i,N/i+i = - S lf S i,N/4+i = S *,i)> the dimer or 
plaquette state has the eigenvalues P = r = 1, while 
the Haldane state P = r = — 1. Therefore, we can dis- 
tinguish these phases by diagonalizing the Hamiltonian 
with twisted boundary in the restricted space specified by 
P = t = ±1. The results for D = 0.4 are shown in Fig. 
|8|. It is found that the two lowest energy levels intersect 
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FIG. 8: Two lowest energies with twisted boundary conditions 
for N = 16 and D = 0.4. The energy of the Haldane state 
(P = t = — 1) and the dimer or the plaquette state (P — 
t = 1) are shown by the open and closed circles. The inset 
shows the Binder parameter (7jv ee ; for J'/J — 0.9. The closed 
and open circles (N = 8 and 16) are calculated with periodic 
boundary conditions. 

each other twice when the ratio of the exchange coupling 
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J'/ J is varied. This implies that the quantum phase 
transitions between the phases with distinct symmetry 
(P = t = ±1) occur twice, being consistent with the 
above predictions based on the VBS analysis. By scaling 
the critical values as (J'/J) C (N) = (J'/J) c (oo) + aN~ 2 , 
we estimate the phase boundaries shown as the crosses 
in Fig. 0, which agree well with those determined from 
the original plaquette chain (N — 16) shown as open 
circles. In contrast, when D is negative, the antifer- 
romagnetic correlation is enhanced, and the systems is 
driven to the Neel ordered phase. To characterize this 
transition, we check the behavior of the Binder param- 
eter, Ujs/ ee i = 1 — 3 ^Q2>2 1 where O is the order pa- 
rameter defined as O = XX - -Sf •. In the inset of 
Fig. ||, the Binder parameters for N = 8 and 16 are 
shown as a function of D with J'/J = 0.9. Since the 
Binder parameter stays invariant with the change of iV at 
the transition point, we can determine the critical value 
(D/J) c ~ —0.07. Consequently, we end up with the 
phase diagram for the spin S = 1 chain with anisotropy 
as shown in Fig. [7[ 
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FIG. 9: VBS picture of the phases for the S = 3/2 model. 

In summary, we have studied the spin-S" quantum spin 
chain with the orthogonal-dimer structure. The obtained 
phase diagram has a rich structure with various type of 
spin-gap phases which are separated by first-order phase 
transitions inherent in fully frustrated systems. Although 
we have given detailed accounts for the S = 1 case, 
it is straightforward to generalize the discussions to an 
arbitrary-spin case. For example, four distinct spin-gap 
phases for S = 3/2 in Fig. || are completely classified 
by the VBS states shown in Fig. ^. Finally we empha- 
size that the first-order phase transitions in the present 
model are triggered by strong frustration, which may cap- 
ture most of essential properties common to this class of 
fully frustrated systems. In particular, we think that 



frustration-induced spin-gap phases obtained here may 
play a key role to clarify the phase diagram of the 2D 
spin-S 1 orthogonal-dimer model. 
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